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Abstract 



Results are reported from a search for non-standard-model Higgs boson decays to 
pairs of new light bosons, each of which decays into the ji^ji^ final state. The 
new bosons may be produced either promptly or via a decay chain. The data set 
corresponds to an integrated luminosity of 5.3 fb^^ of proton-proton collisions at 
\/s = 7TeV, recorded by the CMS experiment at the LHC. Such Higgs boson decays 
are predicted in several scenarios of new physics, including supersymmetric models 
with extended Higgs sectors or hidden valleys. Thus, the results of the search are rel- 
evant for establishing whether the new particle observed in Higgs boson searches at 
the LHC has the properties expected for a standard model Higgs boson. No excess of 
events is observed with respect to the yields expected from standard model processes. 
A model-independent upper limit of 0.78 ± 0.05 fb on the product of the cross section 
times branching fraction times acceptance is obtained. The results are applicable to a 
broad spectrum of models and are also compared with the predictions of two bench- 
mark scenarios as functions of a Higgs boson mass within the range 86-150 GeV/c^ 
and of a new light boson mass within the range 0.25-3.55 GeV/c^. 
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The observation of a new particle lHJIU with a mass near 125GeV/c^ in searches for the stan- 
dard model (SM) Higgs boson [3H3 at the Large Hadron Collider (LHC) raises the critical 
question of whether the new particle is in fact the SM Higgs boson. The precision of the com- 
parisons of the new particle's production and decay properties with the final states predicted 
by the SM will improve with additional data. However, distinguishing a true SM Higgs boson 
from a non-SM Higgs bosons with couplings moderately different from the SM values will re- 
main a challenge. Searches for non-SM Higgs boson production and decay modes are therefore 
particularly timely as they provide a complementary path, which in many cases can allow a 
discovery or rule out broad ranges of new physics scenarios with existing data. 

This Letter presents a search for the production of a non-SM Higgs boson (h) decaying into a 
pair of new light bosons (a) of the same mass, which subsequently decay to pairs of oppositely 
charged muons (dimuons) isolated from the rest of the event activity (h — > 2a + X — > 4/^ + X, 
where X denotes possible additional particles from cascade decays of a Higgs boson). This 
sequence of decays is predicted in several classes of models beyond the SM. One example is 
the next-to-minimal supersymmetric standard model (NMSSM) I|6l- [14| , which extends of the 
minimal supersymmetric standard model (MSSM) [i15h17 | by an additional field (gauge singlet 
under new U{1)pq symmetry) in the Higgs sector of the superpotential. Compared to the 
MSSM, the NMSSM resolves the so-called ^/-problem (generation of the mass parameter }i in 
the Higgs superpotential at the electro weak scale |fT8l ) and significantly reduces the amount of 
fine tuning required lll9ti2T) . 

In the NMSSM, the CP-even Higgs bosons hj and h2 (one of them is the SM-like Higgs boson) 
can decay via hi 2 — > 2ai, where ai is a new CP-odd light Higgs boson p2l - l26 l. The Higgs 
boson production cross section may differ substantially from that of the SM, depending on 
the parameters of a specific model: e.g., a high singlet fraction in hi (or h2) suppresses its 
couplings to SM particles. The new light boson ai couples weakly to SM particles, with the 
coupling to fermions proportional to the fermion mass, and can have a substantial branching 
fraction B{ai V-^V-^) if its mass is within the range 2m^ < nia^ < Inir Il27ll28| . 

Pair production of light bosons can also occur in supersymmetric models with additional hid- 
den (or dark) valleys ||29l - l3T| , which can be motivated by the excesses in positron spectra ob- 
served by satellite experiments 11321 1331 . These dark-SUSY models predict cold dark matter 
with a mass scale of ~ 1 TeV/c^, which can provide the right amount of relic density due to the 
Sommerfeld enhancement in the annihilation cross section arising from a new U{1)d symme- 
try |l34l |35| . In these models, U{1)d is broken, giving rise to light but massive dark photons 
7d that weakly couple to the SM particles via a small kinetic mixing [36 - 38] with photons. The 
lightest neutralino ni in the visible (as opposed to hidden) part of the SUSY spectrum is no longer 
stable and can decay via e.g. ni — )■ 70 + no, where no is a light dark fermion (dark neutralino) 
that escapes detection. The SM-like Higgs boson can decay via h — > 2ni, if > Inin^. The 
branching fraction S (h — )■ 2ni ) can vary from very small to large (bounded by the LHC mea- 
surements of the h — ?► ZZ decay branching fraction). The lack of an anti-proton excess in the 
measurements of the cosmic ray spectrum constrains the mass of 70 to be < GeV/c^ II39I . 
Assuming that 7d can only decay to SM particles, the branching fraction ;B(7d — > H^]^^) can 
be as large as 45%, depending on m-y^ [31J. The Higgs boson production cross section may 
or may not be enhanced compared to the SM, depending on the specific parameters of the 
model. The search described in this Letter was designed to be independent of the details of 
specific models, and the results can be interpreted in the context of other models predicting the 
production of the same final states. 

Previous searches for the pair production of new light bosons decaying into dimuons were 
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performed at the Tevatron with a 4.2 fb data sample ||40| arid more recently at the LHC with 
a 35 pb^^ data sample [41j. Associated production of the light CP-odd scalar bosons has been 
searched for at e+e^ colliders |Il2l|43l and the Tevatron [44J . Direct production of the ai has been 
studied at the LHC [45 1 , but in the framework of NMSSM the sensitivity of these searches is 
limited by the typically very weak coupling of the ai to SM particles. The most stringent limits 
on NMSSM are provided by the WMAP data [46] and LEP searches ll47M49l (rnhj > 86 GeV/c^). 
In the framework of dark SUSY, experimental searches for 70 have focused on the production 
of dark photons at the end of SUSY cascades at the Tevatron [50-52J and the LHC [41 J. Fur- 
thermore, if the newly observed particle at the LHC lU |2l is indeed a Higgs boson, the studies 
of its SM decays will provide additional constraints on the allowed branching fractions for the 
non-SM decays. 

The analysis presented in this Letter uses experimental data collected by the Compact Muon 
Solenoid (CMS) experiment at the LHC. The central feature of the CMS apparatus is a supercon- 
ducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the super- 
conducting solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electro- 
magnetic calorimeter, and a brass/ scintillator hadron calorimeter. The inner tracker measures 
charged particles within the pseudorapidity range |//| < 2.5, where t] = — ln[tan(0/2)] and 
9 is the polar angle with respect to the direction of the counterclockwise proton beam that is 
the z-axis of the CMS reference frame. The tracker provides an impact parameter resolution 
of ~15fim and a transverse momentum (px) resolution of about 1.5% for lOOGeV/c particles. 
Muons are measured in gas-ionization detectors embedded in the steel return yoke. The detec- 
tors are made using the following technologies: drift tubes (\f]\ < 1.2), cathode strip chambers 
(0.9 < \rj\ < 2.4), and resistive-plate chambers (|j/| < 1.6). Matching the muons to the tracks 
measured in the silicon tracker results in a transverse momentum resolution between land 5% 
for px values up to 1 TeV/c. A more detailed description can be found in Ref. [53J. 

The search is performed as a "blind" analysis (the signal region was not used to define the re- 
construction and selection procedures) based on a data sample corresponding to an integrated 
luminosity of 5.3 fb^^ collected with a trigger selecting events containing at least two muons, 
one with px > 17GeV/c and one with px > 8GeV/c. In the offline analysis, events are selected 
by requiring at least one primary vertex reconstructed with at least four tracks and with its z 
coordinate within 24 cm of the nominal collision point. Offline muon candidates are built using 
tracks reconstructed in the inner tracker matched to track segments in the muon system, using 
an arbitration algorithm [54J. The candidates are further required to have at least eight hits 
in the tracker, with the ;\;^/Ndof < 4 for the track fit in the inner tracker (where Ndof is the 
number of degrees of freedom), and at least two matched segments in the muon system. The 
data are further selected by requiring at least four offline muon candidates with px > 8GeV/c 
and \tj\ < 2.4; at least one of the candidates must have px > 17GeV/c and be reconstructed in 
the central region, \rj\ < 0.9. Application of the selection requirements described above yields 
1, 745 events in the data. The trigger efficiency for the selected events is high (96-97%) and 
is nearly independent of the px and t] of any of the four muons. The < 0.9 requirement 
is tighter than that imposed by the trigger, but eliminates significant model dependence at- 
tributable to the reduced trigger performance in the forward region in the presence of multiple 
spatially close muons. This rj requirement causes an overall reduction in the analysis accep- 
tance of about 20%, as obtained in a simulation study with one of the NMSSM benchmark 
samples used in the analysis. 

Next, oppositely charged muons are grouped into dimuons (a muon may be shared between 
several dimuons) if their pairwise invariant mass satisfies m^/^, < 5GeV/c^ and if either the fit 
of the two muon tracks for a common vertex has a fit probability greater than 1% or the two 
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muon tracks satisfy the cone size requirement AK(f/+, }i ) = \J — f]^- )^ + (^^+ — (ppi- )^ < 
0.01, where (p is the azimuthal angle in radians. The AR requirement compensates for the re- 
duced efficiency of the vertex probability requirement for dimuons with very low mass (m^^ > 
2m^), in which the two muon tracks are nearly parallel to each other at the point of closest 
approach. 

Once all dimuons are constructed, only events with exactly two dimuons not sharing common 
muons are selected for further analysis. There is no restriction on the number of ungrouped 
(orphan) muons. Assuming that each dimuon is a decay product of a new light boson, we 
require that the two dimuons have invariant masses in the range 0.25-3.55 GeV/c^. We recon- 
struct Zpipi, the z coordinate of the dimuon system at the point of the closest approach to the 
beam line, using the dimuon momentum measured at the common vertex and the vertex po- 
sition. We ensure that the two dimuons originate from the same pp interaction by requiring 
\zptpt^ — Zjiji2 1 < 1 mm. This selection yields 139 events in data and it is fully efficient for signal 
events while reducing the probability of selecting rare events with dimuons from two separate 
primary interactions. 

To suppress SM background, we require that the dimuons be isolated from other activity in the 
event, using the criterion Jgum < 3 GeV/ c, where the isolation parameter of the dimuon system 
hum is defined as the scalar sum of the transverse momenta of all additional charged tracks 
with px > 0-5 GeV/c within a cone of size AR = 0.4 centered on the momentum vector of the 
dimuon system. Tracks used in the calculation of hum must also have a z coordinate at the point 
of the closest approach to the beam line that lies within 1 mm of the z coordinate of the dimuon 
system. The hum selection yields three events in data and it suppresses the contamination from 
bb production by about a factor of 40 (measured in data) while rejecting less than 10% of the 
signal events (obtained from the simulation study). 

Finally, we require that the invariant masses of the two reconstructed dimuons are compatible 
with each other within the detector resolution \mi — mz] < 0.13 GeV/c^ + 0.065 x (nii + m2)/2, 
where mi = m^^j and mi = mj,^^- The numerical parameters in this last requirement cor- 
respond to at least five times the size of the core resolution in dimuon mass, including the 
differences in resolution in the central and forward regions. The inefficiency of this W2i ~ W22 
selection is less than 5% per event; it is due to QED final-state radiation and is unrelated to the 
detector resolution. 

To illustrate the performance of the analysis to select possible signal, we use the two benchmark 
models introduced earlier. The NMSSM samples are simulated with the pythia 6.4.26 event 
generator [55 1 using MSSM Higgs boson production via gluon-gluon fusion gg — > H^ggj^, 
where the Higgs bosons are forced to decay via H^ggj^ —> 2A^c,g^- The masses of H^ggj^ and 
^MSSM to the desired values for the hi mass (within the range 86-150 GeV/c^) and ai 

mass (within the range 0.25-3.55 GeV/c^), respectively. Both A^ggj^ bosons are forced to decay 
to a pair of muons. The dark-SUSY samples are simulated with the MadGraph 4.5.2 event 
generator Il56l using SM Higgs boson production via gluon-gluon fusion gg — > hgM, where 
the mass of hg^ is set to the desired value for the h mass (within the range 90-150 GeV/c^). 
The BRIDGE software [57] was used to implement the new physics model that forces the Higgs 
bosons hgM to undergo a non-SM decay to a pair of neutralinos ni, each of which decays ni — > 
riD + 7d/ where mm = lOGeV/c-^, ntn^ = 1 GeV/c^ and m-y^ = 0.4GeV/c^. Both dark photons 
7d are forced to decay to two muons, while both dark neutralinos n^ escape detection. In each 
case the narrow width approximation is imposed by setting the widths of the Higgs bosons 
and dark photons to a small value (lO^"' GeV/c^). 
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Table 1: Event selection efficiencies £{m^^,m^^) with statistical uncertainties for the NMSSM 
benchmark scenario, as obtained from simulation. 
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Table 2: Event selection efficiencies £{m-^,m^^) with statistical uncertainties for a dark-SUSY 
benchmark model, as obtained from simulation. 
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All events of the benchmark signal samples are processed through a detailed simulation of 
the CMS detector based on Geant4 [58J and are reconstructed with the same algorithms used 
for data analysis. Tables [l] and |2] show the event selection efficiencies efuu obtained using the 
simulated signal events for these two benchmark scenarios. To provide a simple recipe for 
future reinterpretations of the results in the context of other models, we separately determine 
*gen/ the geometric and kinematic acceptance of this analysis calculated using generator level 
information only. It is defined with the criteria that an event contains at least four muons 
with px > 8GeV/c and \r]\ < 2.4, with at least one of these muons having pj > 17GeV/c 
and I?/ 1 < 0.9. Tables [l] and |2] also show ctgen along with the ratio ^fuu/^i^gen- The systematic 
uncertainty for the ratio efuu/is^gen is 7.4%, for agen is 3.0% and for efuu is 8.0% as described 
below. The model independence of the ratio permits an estimate of the full event selection 
efficiency of this analysis for an arbitrary model predicting the same signature with two pairs 
of muons. The acceptance agen is calculated using a suitable event generator. The full efficiency 
efuu is then calculated by multiplying agen by the ratio £iua/ *gen = 0.74 ± 0.05. 

The background contributions after final selections include bb and direct ]/xp pair production 
events. The leading part of the bb contribution is due to b-quark decays to pairs of muons via 
double semileptonic decays or resonances, i.e. co, p, (p, ]/xp. A smaller contribution comes from 
events with one real dimuon and a second dimuon with a muon from a semileptonic b-quark 
decay and a charged hadron misidentified as another muon, due to the incorrect association of 
the track of the charged hadron with the track segments from a real muon in the muon system. 
The contribution of other SM processes has been found to be negligible, e.g. the low mass 
Drell-Yan production is heavily suppressed by the requirement of additional muons. 

The direct ]/xp pair production contribution is estimated using the simulation scaled to data in 
a few regions of the invariant mass of the j/xp pair. The comparison yields an estimate for the 
number of direct ]/tp pair produced events satisfying all analysis criteria of 0.3 ± 0.3 events. 

The bb background contribution is modeled as a two dimensional (2D) template B^^{m-[, 1172) in 
the plane of the invariant masses of the two dimuons in the selected events, where nii always 
refers to the dimuon containing a muon with pj > 17GeV/c and \f]\ < 0.9. For events with 
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nrii (i = 1 ,2) [GeV/c^] [GeV/c^] 

Figure 1: Left: Comparison of the data (solid circles) failing the mi ~ m2 requirement in the 
control sample (no isolation requirement is applied to reconstructed dimuons) with the predic- 
tion of the background shape model (solid line) scaled to the number of entries in the data. The 
insets show the B17+8 and Bs+s templates (solid lines) for dimuons obtained with background- 
enriched data samples. Right: Distribution of the invariant masses m-[ vs. m2 for the isolated 
dimuon systems for the three events in the data (shown as empty circles) surviving all selec- 
tions except the requirement that these two masses fall into the diagonal signal region mj ~ m2 
(outlined with dashed lines). The background expectation (as indicated by the intensity (color 
online) of the shading) is a sum of the bb and the direct J/1/7 pair production contributions. 

both dimuons containing such a muon, the assignment of nii and m2 is random. As each b 
quark fragments independently, we construct the template describing the 2D probability den- 
sity function as a Cartesian product Bi7+8("^i) x B8+8(?"2)/ where the B17+8 and B8+8 templates 
model the invariant-mass distributions for dimuons with or without the requirement that the 
dimuon contains at least one muon satisfying pj > 17GeV/c and \ri\ < 0.9. This distinction 
is necessary as the shape of the dimuon invariant mass distribution depends on the trans- 
verse momentum thresholds used to select muons and whether the muons are in the central 
(I //I < 0.9) or in the forward (0.9 < < 2.4) regions, owing to the differences in momentum 
resolution of the barrel and endcap regions of the tracker. The Biy+s shape is measured using 
a data sample enriched in bb events with exactly one dimuon and one orphan muon under the 
assumption that one of the b quarks decays to a dimuon containing at least one muon with 
Pt > 17GeV/c and \ f]\ < 0.9, while the other b quark decays semileptonically resulting in an 
orphan muon with pj > 8 GeV/ c. For the B8+8 shape, we use a similar sample and procedure 
but only require the dimuon to have both muons with pj > 8 GeV/c, while the orphan muon 
has to have px > 17GeV/c and \f]\ < 0.9. Both data samples used to measure background 
shapes are collected with the same trigger and with kinematic properties similar to those bb 
events passing the selections of the main analysis. These event samples do not overlap the 
sample containing two dimuons that is used for the main analysis, and they have negligible 
contributions from non-bb backgrounds. The B17+8 and B8+8 distributions, fitted with a para- 
metric analytical function using a combination of Bernstein polynomials I.59J and Crystal Ball 
functions l|60ll describing resonances, are shown as insets in Fig. [T] (left). Once the B^^[mi,m2) 
template is constructed, it is used to provide a proper description of the bb background shape 
in the main analysis. 
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To validate the constructed B^^{mi,m2) template, we compare its shape with the distribution 
of the invariant masses mi vs. m2 from events obtained with all standard selections except 
the requirement that each of the two reconstructed dimuons is isolated. Omitting the isolation 
requirement provides a high-statistics control sample of events with two dimuons highly en- 
riched with bb events. To avoid unblinding the search, the diagonal signal region is excluded 
in both the data and the template, i.e. the comparison has been limited to the data events that 
satisfy all analysis selections but fail the nii ~ m2 requirement. Distributions of nii and m2 are 
consistent with the projections of the B^^{m-[,m2) template on the respective axes normalized 
to the number of events in the data control sample. The sum of the mi and m2 distributions 
agrees well with the sum of the template projections as shown in Fig.[l](left). 

Another cross-check has been performed using data events which satisfy all analysis selec- 
tions except that the isolation parameters of each dimuon system have been required to satisfy 
3 GeV/ c < /sum < 8 GeV/ c, which removes potential signal events (the signal selections require 
hum < 3GeV/c for each dimuon). These selections yield four events in the off-diagonal side- 
band region of (mi, m2) plane. Normalizing the background distribution to these four observed 
events, we predict 0.9 ± 0.4 bb events in the diagonal region, consistent with no events being 
observed there. 

To normalize the constructed B^^{mi, m2) template, we use the data events that satisfy all anal- 
ysis selections except for the mi ~ m2 requirement. These selections yield three events in the 
off-diagonal sideband region of (mi,m2) plane, which corresponds to an expected number of 
0.7 ± 0.4 bb events in the diagonal signal region. This result is shown in Fig. [l] (right) where 
the three events in the off-diagonal sidebands of the (mi, m2) plane are shown as empty circles. 
The distribution of the total background expectation in the (mi, m2) plane (a sum of the bb and 
the direct ]/^p pair production contributions) is also shown by the intensity of the shading in 
Fig. [1] (right). 

Selection efficiencies of muon and dimuon reconstruction, trigger, and isolation criteria are 
obtained with the simulation and have been corrected by using the scale factor of 0.91 (per 
event), obtained by comparing data and simulation. Systematic uncertainties include the un- 
certainty in the LHC integrated luminosity of the data sample (2.2%) [61 1, uncertainties related 
to parton distribution functions (PDF) and the knowledge of the strong coupling constant as 
(3%), offline muon reconstruction (5.7%), dimuon reconstruction that accounts for effects re- 
lated to overlaps of muon trajectories in the tracker and in the muon system (3.5%), trigger 
(1.5%), and the dimuon mass shape, which affects the efficiency of the requirement that the 
two dimuon masses are compatible (1.5%). All uncertainties quoted are relative to the final 
analysis selection efficiency per signal event. The PDF uncertainties are estimated by compar- 
ing the CTEQ6.6 [62] with NNPDF2.0 [|63| and MSTW2008 [64J sets following the PDF4LHC 
recommendations Il65l . Varying the QCD renormalization/ factorization scales has been found 
to have a negligible effect. The total systematic uncertainty in the selection efficiency is 8.0%. 

When the final data was unblinded, no events were observed in the signal diagonal region, as 
illustrated in Fig. [l] (right). The expected number of background events in the signal diagonal 
region is 1.0 ± 0.5. This number includes contributions from the bb and the direct J/i/? pair 
production. 

We interpret these results in the context of the dark-SUSY and the NMSSM scenarios and obtain 
95% confidence level (CL) upper limits, using a Bayesian prescription. For the dark-SUSY 
model, the limit is set on (7(pp — > h — > 2ni — > 27d +2nD) x B^{'Yd — ^ 2f/) as a function 
of mil. This limit is shown in Fig. |2] (left) for m-y^ = 0.4GeV/c^, which yields the branching 
fraction B{jd — ^ 2^) close to its maximum, reaching approximately 45% [31]. In the case of the 
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Figure 2: Left: The 95% CL upper limit as a function of mh, for the dark-SUSY scenario, on 
(7"(pp — > h — > 27d +2nD) x B^{jd — ^ 2f/) with niy^ = 0.4GeV/c^ (dashed curve). As an 
illustration, the limit is compared to the rate obtained using a toy model with SM Higgs boson 
production cross section ^(pp — > h) = (Tsm, B{h — > 2ni) = 1%, B{ni 70 H-no) = 50%, 
and i3(7D —> 2}i) = 45% (solid curve). The chosen i5(7D —> Iji) is taken from ||3T]| . Right: 
The 95% CL upper limits as functions of m^^, for the NMSSM case, on c^(pp hi 2 2ai) x 
S^(ai — > with = 0.25GeV/c^ (dashed curve), = IGeV/c^ (dash-dotted curve) 
and maj = 3.55GeV/c^ (dotted curve), compared to the rate obtained using a toy model with 
C7"(pp hi) = cTgM, S(hi 2ai) = 3%, and jB(ai — > 2n) = 7.7% (solid curve). The chosen 
6(ai ^ 2f/) is taken from 1281 . 

NMSSM, the 95% CL upper limit is set for a (pp hi,2 ^ 2ai) x S^(ai ^ 2^) as a function of 
for three choices of maj as shown in Fig.|2](right) and as a function of m^-^ for three choices of 
Mthj as shown in Fig. [s] (left). We also set the 95% CL upper limit for B (hi^2 — > 2ai) x 5^(ai 
2}i) as a function of maj for three choices of assuming (7(pp — > hi) = cr(pp h2) = 
(7gM(125GeV/c^) as shown in Fig. [s] (right). As m^^ is unrestricted for any given m^^, we use 
^fuii("^h2) = ^fuii(^hi) to simplify the interpretation. This is conservative since €{ua_{m^^) > 
efuii(mhi) if > mhi, for any nia^. 

In the representative scenarios, for any fixed combinations of fWi, and both the Higgs bo- 
son production cross section and the branching fractions can vary significantly, depending on 
the model parameters. In the absence of broadly accepted "model slopes" for the NMSSM 
or the dark-SUSY scenarios, we compare the obtained experimental limits with a few bench- 
mark "toy" models, assuming the SM Higgs boson production cross section as a function of the 
Higgs boson mass ||66J. We further use S(7d —> 2ji) = 45% [31 1 for the dark-SUSY model and 
we assume B{ni — > 7d + no) = 50%, allowing for other possible decays. We use B{ai 2}i) 
as a function of JWaj, obtained using the NMSSM calculation (hadronization effects in region of 
low JMaj are not included) in [28 1 for tanjS = 20. The branching fraction ;B(ai — > 2f/) is equal 
to 7.7% for mai ~ 2GeV/c2. Finally, we choose B{h 2ni) = 1% and B{hi^2 2ai) = 3%, 
which yields the toy model predictions for rates of the dimuon pair events comparable to the 
experimental limits. With these parameter values, the branching fractions include the Higgs 
boson couplings to the new light bosons as well as the reduced couplings of Higgs bosons to 
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Figure 3: Left: The 95% CL upper limits as functions of nia^, for the NMSSM case, on cr(pp — > 
hi,2 ^ 2ai) X i3^(ai -> 2}i) with = 86GeV/c2 (dashed curve), rnhj = 125GeV/c^ (dash- 
dotted curve) and m^^ = 150GeV/c^ (dotted curve). The limits are compared to the rate (solid 
curve) obtained using a toy model with cr(pp hi) = cr(pp — > h2) = c7sM(125GeV/c^), 
^{^1,2 ^ 2ai) = 3%, and S(ai ^ 2]d) as a function of m^^ which is taken from ||2H1 for NMSSM 
parameter tan ^ = 20. Right: The 95% CL upper limits on B{hi^2 ^ 2ai) x B^{ai 2}l) with 
njhi = 86 GeV/c^ (dashed curve), m-^^ = 125 GeV/c^ (dash-dotted curve) and /Whi = 150 GeV/c^ 
(dotted curve) assuming cr(pp — > hi) = cr(pp h2) = crg]y[(125GeV/c^). The limits are 
compared to the branching fraction (solid line) obtained using B (hi^2 — > 2ai ) = 3% and B (ai — > 
2}i) as a function of nia^ which is taken from p8|. 

the SM particles in these models, which affect the Higgs boson production cross section. The 
sensitivity of this search can be compared to that of a similar analysis performed at the Teva- 
tron [40 1 after rescaling with the ratio of the Higgs boson cross sections at the LHC and the 
Tevatron. If plotted in Fig. |2] and Fig. |3] (left), the Tevatron results would have excluded rates 
above ~ 130 fb, an order of magnitude in sensitivity less. 

For an arbitrary non-SM scenario predicting the signature investigated in this Letter, the results 
can be presented as the 95% CL limit cr(pp ^ 2a + X) x B'^{di -> 2/^) x ttgen < 0.78 ± 0.05 fb, 
where agen is the generator level kinematic and geometric acceptance described earlier. The 
calculation uses integrated luminosity C = 5.3 fb^^ and takes the ratio efuii/ *gen = 0.74 ± 0.05, 
which includes the systematic uncertainties and covers the variation in the ratio over all of the 
benchmark points used (see Tables [l] and |2|. The limit is not applicable to models where the 
new light bosons are typically non-isolated or have substantial lifetime. The efficiency of the 
selections in this analysis abruptly deteriorates if the light boson's decay vertex is more than 
~ 4 cm away from the beamline. 

In summary, no excess is observed in the data with respect to the SM predictions. We find no 
evidence of non-SM decay modes of a Higgs boson into pairs of new light bosons of the same 
mass, which subsequently decay to pairs of oppositely charged muons (h ^ 2a + X ^ 4/^ + X) 
for Higgs boson masses in the range 86 < < 150 GeV/c^ and for the new light-boson masses 
in the range 0.25 < wJa < 3.55GeV/c^ using data collected by the CMS experiment in proton- 
proton coUisioris at ^/s = 7TeV, corresponding to an integrated luminosity of 5.3 fb^^. The 
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analysis has been designed as a quasi-model-independent search allowing interpretation of its 
results in the context of a broad range of non-SM scenarios predicting this type of signature. 
In the context of the NMSSM and one of the SUSY models with hidden valleys this search 
provides the best experimental limits to date, significantly surpassing the sensitivity of earlier 
similar searches performed at the Tevatron. 
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